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THE WEIL-PETERSSON HESSIAN OF LENGTH ON 
TEICHMULLER SPACE 



MICHAEL WOLF 

Abstract. We present a brief but nearly self-contained proof of a for- 
mula for the Weil-Petersson Hessian of the geodesic length of a closed 
curve (either simple or not simple) on a hyperbolic surface. The formula 
is the sum of the integrals of two naturally defined positive functions over 
the geodesic, proving convexity of this functional over Teichmuller space 
(due to Wolpert (1987)). We then estimate this Hessian from below in 
terms of local quantities and distance along the geodesic. The formula 
extends to proper arcs on punctured hyperbolic surfaces, and the es- 
timate to laminations. Wolpert's result that the Thurston metric is a 
multiple of the Weil-Petersson metric directly follows on taking a limit 
of the formula over an appropriate sequence of curves. We give further 
applications to upper bounds of the Hessian, especially near pinching 
loci, and recover through a geometric argument Wolpert's result on the 
convexity of length to the half-power. 



1. Introduction 

One of the foundations of modern Teichmuller theory is Wolpert's [W6187] 
theorem that the function on Teichmuller space that records the geodesic 
length of a simple closed curve is convex with respect to the Weil-Petersson 
metric. That paper provided a lower bound for the Hessian of the length 
function; our purpose here is to present a brief derivation of a concise formula 
for the Hessian in terms of natural objects on the surface associated to the 
curve and the tangent vectors to Teichmuller space. 

To explain this result and add some context, we fix some terminology and 
notation. Let S be a smooth closed surface of genus g, and let T(S) be the 
Teichmuller space of (isotopy classes of) marked hyperbolic structures on 
S. Let [7] be a free homotopy class of closed curves, not necessarily simple, 
a slight generalization of the setting in which Wolpert worked. Typically 7 
will denote the representative of [7] that is geodesic with respect to a given 
metric g, with the context making it clear whether 7 is the immersed curve 
on the surface or an immersion from a circle into the surface. 

For each hyperbolic surface (S,g), there is a geodesic representative 7 = 
7 g of [7]. By the uniformization theorem, as each point in T{S) is rep- 
resented by a unique hyperbolic metric, say g, the length of the geodesic 
representative 7 of [7] defines a function t = £ 7 = ^ 7 ([g]) on T(S). We 
investigate the second derivative of that function. 



2000 Mathematics Subject Classification. Primary 30F60. 

Partially supported by NSF grants DMS-0139877 and DMS-0505603. 

1 



2 



M. WOLF 



The Hessian of a function is well-defined once there is a background met- 
ric. There are many metrics on T(S); among the more basic is the Weil- 
Petersson metric. Representing the tangent space to Teichmiiller space as 
the space of Beltrami differentials which are harmonic with respect to the 
hyperbolic metric g representing a point in T(S), the Weil-Petersson metric 
is the L 2 metric on that space with respect to the hyperbolic area form dA g . 

The first goal of this paper is to write a formula for the Weil-Petersson 
Hessian of the function. 

1.1. Statement of the Formula. In this subsection, we describe our for- 
mula in the simplest case; we discuss the situation of the second derivative 
of the length function I of a closed curve, with the derivative taken along a 
Weil-Petersson geodesic T. Even in this case, we require some notation. 

Let T = T(t) a Weil-Petersson geodesic arc; the class [7] is represented 
by the r(t)-geodesic j t . The tangent vector to Teichmiiller space at T(0) is 
given by a harmonic Beltrami differential, say fi = ^. 

We can extend r(0)-Fermi coordinates along the curve 70 to complex 
coordinates in a neighborhood of (a subarc of) 70. In terms of those co- 
ordinates, the quantity — = Im^ is well-defined. Let denote the 
solution to the ordinary differential equation 



where here the geodesic is represented by a vertical line in the Fermi coor- 
dinate patch with a parametrization given by arclength. 

This is enough terminology so that we may state our main result as the 

Theorem 1.1. Along the Weil-Petersson geodesic arc the second vari- 
ation 4pj£ of the length £(t) = L(T(t), [7]) is given by 



1.2. Applications and Extensions. The formula (1.2) lends itself to ap- 
plications. Using well-known techniques for solving and estimating solutions 
of the relevant differential equations, we obtain results in a number of dif- 
ferent directions. 

Of course, it is immediately apparent that this Hessian is positive definite, 
even for curves which are not simple. The first summand is written in terms 
of 2(A — 2)~ 1 ^i-, a ubiquitous term in Teichmiiller theory whose definition 



(1.1) 




Im<3? 



(1.2) 
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involves the global geometry of (S,g). In Lemma 5.1, we find a locally 
defined lower bound for this quantity. Thus, when combined with the second 
expression for the second term, we obtain an estimate defined only in terms 
of local quantities or distance along the geodesic. 

In addition to providing a means to estimate the Hessian from below, we 
also obtain an improvement on the basic convexity result: from an easily 
derived formula for the gradient of length, we will see (Corollary 8.1) im- 

2 

mediately that l'* is convex on all of Teichmuller space. A recent theorem 
[Wol08] of Wolpert is that this may be improved to £2 being convex on 
all of Teichmuller space; we provide a proof for that as well. The proof 
is geometric in the sense that it hinges on a comparison of two harmonic 
diffeomorphisms of a cylinder. 

It is straightforward to extend our derivations both to laminations and 
to proper geodesic arcs which connect cusps on a hyperbolic cusped surface. 
The explicit nature of formula (1.2) also lends itself to estimates from above, 
leading to estimates on the Weil-Petersson connection near the pinching 
locus. These were also recently obtained (and announced some time ago) 
by Wolpert [W0IO8]. 

Our final application is a new proof of Wolpert 's [Wol86b] proof that 
the Thurston metric is (a multiple of) the Weil-Petersson metric. We con- 
sider a sequence {j n } of curves whose geodesic representatives are becoming 
equidistributed in the unit tangent bundle T 1 (S', g). Thurston observed that 
linin^oo Hess^ 7n would be a positive definite quadratic form on T^ S ^T(S): 
by taking a limit of the right-hand side of formula (1.2), we see that this 
Thurston metric is a multiple of the Weil-Petersson metric. Wolpert's argu- 
ment followed a more quasiconformal analytic tradition, while this derivation 
is more Riemannian in perspective. 

1.3. Organization of the paper. We organize the paper as follows. In 
Part 1, we derive Theorem 1.1. Computations in Teichmuller theory often re- 
quire fixing a gauge; here we find it convenient to vary hyperbolic structures 
(S,gt) under the condition that the identity mapping id : (S,go) — ► (S,gt) 
is harmonic. This requirement gives the prescribed curvature equation a 
particularly convenient form, and our derivation begins with a sketch of a 
useful computation from [Wol89]. The rest of the derivation is self-contained, 
occupying sections 2-4. 

Part 2 of the paper contains the applications and extensions. In section 
5, we find a lower bound for the integrand in the first term; this leads to 
an estimate (Corollary 5.2) for the Hessian in terms of an integral along the 
curve of local quantities. Section 6 is devoted to a derivation of our results 
in the setting of a curve which connects ends of a (complete) hyperbolic 
punctured surface. Here, while the length of an arc is infinite, the Hessian 
of a regularized version of the length is positive and finite; section 7 extends 
our work from curves to laminations. In section 8, we use the formulae to 
give some geometric estimates: we quickly prove that not only is the length 
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£ of curves convex, but so is l'i . We then give a longer geometric argument 
that £2 is also convex. We give a general upper bound for the Hessian, 
as well as estimates for the Weil-Petersson connection near the Deligne- 
Mumford compactification divisor. Finally, in section 9, we take a limit of 
formula (1.2) in Theorem 1.1 over curves that are becoming equidistributed 
to recover the result that the Thurston metric is a multiple of the Weil- 
Petersson metric. 

Acknowledgments. The author is grateful to David Dumas for conversa- 
tions, especially on lengths of laminations; Yair Minsky for pointing out the 
question on convexity of arcs and for encouragement to find cleaner formu- 
lae; Scott Wolpert for comments, particularly on section 8.1; and especially 
Zheng (Zeno) Huang for very careful reading, criticism, the location of some 
mistakes, and suggestions for improvements. 

Part 1. A formula for the Weil-Petersson Hessian of length 

We begin with a brief background discussion of the some the theory of 
Teichmuller space and the Weil-Petersson metric we will need. Tangent 
vectors to Teichmuller space at a point (S, g) are represented by 'harmonic 
Beltrami differentials' of the form /i = *, where g is a hyperbolic metric on S 
and <E> is a holomorphic quadratic differential on (S,g). The Weil-Petersson 
inner product of two such tangent vectors is the L 2 inner product 

(1.3) (-,-) = Re / dArea q . 

\9 9 I Js 9 9 

Much is now known about this metric: by means of an introduction to 
the subject, the Weil-Petersson metric is not complete [Chu76] [W6175], it 
is Kahler [AI1I6I], negatively curved with good expressions [Roy] [Tro86] 
[Wol86a] for and estimates [Hua05] [LSY04] of the curvatures, it is quasi- 
isometric to the pants complex [Bro03], the isometry group is exactly the 
(extended) mapping class group [MW02]. Of course, the stimulus for this 
article and an important ingredient in some of the results above is that the 
Weil-Petersson metric is geodesically convex [W6187]. 

2. The second derivative of length in space and time. 

We are interested in computing the second variation of geodesic length of 
a curve along a Weil-Petersson geodesic. We imagine the setting as a fixed 
differentiable surface S equipped with a family of metrics gt, and on this 
surface there is a family of curves ■ The curves are all freely homotopic 
and may or may not be simple. The defining equation is that the curves jt 
are ^-geodesies; we shall shortly write that equation out in coordinates. 

To begin though, we separate the overall second variation of length into 
a term that refers only to the second variation of the metric gt and a term 
that refers only to the second variation of the curve jf This separation is 
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quite standard for a variational functional. Write the length of 7 S in the 
metric gt as L(gt,j s ). Then, in this language, the geodesic equation takes 
the form, for all t, 



(2.1) 



d_ 

ds 



L(gt,is 



s=s 



d_ 
~d~s 



o. 



if 7 S0 is a ^-geodesic, and Jj7 s is an infinitesimal variation of curves through 

The second variation of length of the ^-geodesies is given by 
(2.2) 
d 2 

—L(g ult ) = £>?iL(so,7o)[5,5] + 2D 2 12 L(g , 7o M[i} + D 2 22 L(g 0l 7 o)[7,7]- 

where g = ^g t and 7 = Of course, if jt is a ^-geodesic, then we write 
the geodesic equation (2.1) above in this notation as 



D 2 L(g t , lt )[y] =0 



and so 



d 



= -D 2 L(g t , lt m 
= D 1 D 2 L(g ,j )[g}[ A f] + D 2 D 2 L(g , 70) [7, 7]. 



Thus, 
(2.3) 



(2.4) 



Dl 2 L i90, To) [9] [7] = -D% 2 L(g , 70) [7,7]. 
Substituting (2.3) into (2.2) yields that 
d 2 



dt 2 



L(gt,lt) = D u L(go,7o)[g,g] - D 22 L(g , 7o ) [7] [7]. 



Some remarks on this equation (2.4) are in order. First, note that the 
term D 22 L(go, 70) [7] [7] is non-negative, as the surface (S,go) is negatively 
curved; indeed, this second variation term is positive unless the vector field 
7 is tangent to the curve 70. Thus our task is to prove that the first term 
D 2 1 L(go, 7o)[<7, g] is larger than the second term. 

In the next sections, we evaluate the terms D^L and D 22 L via different 
methods. 



3. Second variation of arclength of 70 in a family of metrics 

We briefly recall the computational scheme of [Wol89]. Let $ G QD(go) 
denote a quadratic differential, holomorphic with respect to a conformal 
metric go . Then we may consider a family of metrics on S decomposed by 
type as 

/ t 2 |$| 2 \ — 

(3.1) g t = t^dz 2 + g \H(t) + ^T^y J dzdz + t$dz 2 . 
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Here z is a conformal coordinate for (S,go). It is straightforward to check 
[SY78] that the metric gt is hyperbolic if 



(3.2) 



A 90 logW«)=2W«)-|ffl!-2. 



(Of course, the pullback of a hyperbolic metric by a diffeomorphism is hy- 
perbolic, and so we might imagine that if we pullback gt by a family of 
diffeomorphisms ipt, then the result ipt9t would also be hyperbolic. Here we 
have chosen a gauge by requiring that the identity map id : (S, go) — > (S, gt) 
is harmonic.) 

We are interested in second variations. Differentiating twice and applying 
the maximum principle to the first derivative (see [W6189] for an expanded 
description) yields 



(3.3) 



n 
n 



dt 
dt 2 



t=o 



H(t) = 

H{t) = -2(A- 2)" 1 



2|$| 



t=o % 

We observe that — 2(A— 2) _1 is a positive operator and so H > 0. Combining 
(3.1) and (3.3), we conclude that 

(3.4) 

9t = g{t) = godzdz + t(®dz 2 + ®dz 2 ) 

2|$| 2 ' 



+ r/2 



2|<1>| 2 
9o 



9l 



g$dzdz + 0(t 



Now use that 



DnL(g ,io)[g,g] = 



dt 2 



L(gt,jo) 



dt 2 



t=Q -> 70 



gt- 



Substituting (3.4) into this last integral with a choice of coordinate so that 
70 is a line {Re z = const} and differentiating under the integral symbol 
then yields 

(3.5) D 2 1 L(g , l0 )[g,g] = [ -I( 5o )- 3 / 2 (2 Re <f) 2 



2|$| 2 /A , ,2\$\ 2 
^-2(A-2)" 1 ^ 



9o 



9l 



fi 



(Im$) 
~9T 



2(A-2) 



-i 



I* 



21 



9o J 



since |<I>| 2 — (Re<3?) 2 = (Im<I>) 2 . Both terms are positive, and so we see that 
this expression is positive, as we expected (and needed if the expression 
D\ V L — D 2 2 L is to be positive). 
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Remark. As an easy model of this method, we quickly reproduce a formula 
for the first variation of length. (We'll have use of this expression in later 
sections.) 

We compute the first derivative of length i along T(t) to be 
(3-6) j t £ 1 (T(t)) = D 1 L(g t , lt )[g}+D 2 L(g t , lt )[j]. 

Of course, as 70 is a geodesic, the second term D2L{g t , 7t)[7] = and 
from (3.1) and (3.2), we find 

< 3 - 7 > sW» = J t j^ 

d n 

m9t 



/ 70 2go 

f Re$ , 
= / as, 

J 70 5o 

by (3.4), concluding the computation. 

4. Second variation of ^o-arclength of the family j t of 

Cft-GEODESICS 

Our next step is to evaluate the term D^Lfao, 70) = ^L(go, 7t)[7,7] in 
(2.4). 

It is of course standard (see for example [Spi79b]) that if V is the varia- 
tional field of a family of curves through a geodesic 70, then 

-,2 



(4.1) = I 



dV 1 
ds 



-K\V\ 2 ds 



where K = K(s) denotes the Gaussian curvature of the surface at the point 
7o(s). To compare this formula (4.1) to (3.5), we will need to find an ex- 
pression for V in terms of the quadratic differential <£. That is the main 
goal of this section. 

Of course, the defining equation of jt is that it is a geodesic, or equiva- 
lently that its geodesic curvature vanishes. We write this schematically, in a 
similar way that we write the length L = L(g t ,~ft), as a function k = K(gt,7t) 
of a metric and a curve: 

(4.2) «Gft,7t) = 0. 

Differentiating in t, we find that 

(4-3) ^(#o,7o)[7] = -^(50, 7o) [3]- 

As expected, the left-hand side of (4.3) is the classical Jacobi operator 

+-^) but the right hand side will involve the first derivatives of gt, 

i.e. the metric go and the quadratic differential Our next task will be to 
find an expression for the solution 7 to (4.3). 
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4.1. The inhomogeneous Jacobi equation. We first expand the right 
hand side of (4.3). We pick conformal (Fermi) coordinates z = x + iy so 
that the geodesic 70 is described by {x = const}. These coordinates are a 
bit unusual in that the geodesic may repeatedly visit the same points on the 
surface: it's possibly better to regard the geodesic as embedded in the unit 
tangent bundle T l M with z = x + iy its projection to the surface S. 

Then, invoking a coordinate expression for n (see [Opr04], consistent with 
definitions in [Spi79a] and [Spi79b]), we have 



d -K(go,jo)[g] = 




1 Vdetff(t) 



dt^> dt 522 (t)3/2 

where g(t) = gijit) is defined by (2.1) as 
/<7ii(t) 512 (t) N 

(4-4) g(t) = 

\922{t) 922 (t) / 

'g + 2tRe$ -2tlm$ 



E F 
F G 



+ 0{t 2 
:st func 

end as it actually simplifies some of our notation; for example, we write 



-2£lm$ 3o-2tRe$, 
We include the very classical notation for the first fundamental form at the 



(4.5) K(t) = -rUt) VE ^ /2 F2 -, 

where here of course the variables E = E(t), F = F(t), and G = G(t) all 
depend on t. Now, in this language, suppressing some of the dependence on 
t, we have 

, , . „i . . *2GFi, — GG T — FGu 

(4.6) r> 2(() = l {m _ n 1 . 

Since k(0) = and F(0) = 0, we find that 

(4.7) ^ 922 (0) = 1-G = on 70 . 

We differentiate (4.6) in t and use (4.7) and F(0) = to find that 
J t Tl 22(t) = ^2{-45o(Im$) y + 2g (Re<S>) x + 2(Im $)( 5o ) y } 

= A{-2fld(Im*)„ + 2(Im*)( go ) 1/ }; 

here the last equality folows from the Cauchy-Riemann equations for the 
real and imaginary parts of the holomorphic quadratic differential <£. We 
conclude that 
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(4.10) V yy - V = — lm/I. 



Combining (4.3), (4.5)and (4.8) yields the equation we will focus on: 

(4.9) £y_y« 

oy l dy { g 

Remark. It is easy to compute that the Beltrami differential tangent to 
our deformation is given by fi = J^. In that language, our equation (4.9) 
becomes 

d_ 

dy 

4.2. The Primitive of the Variation Field. The right-hand side of (4.3) 
is the derivative of the basic quantity — ^sl* appearing in (3.5). This term 

provides a link between the two different terms D^L and D 22 L of the basic 
expression (2.4) for the Hessian of length. To find the final formula for the 
Hessian of the length function, we consider the primitive of V = 7 along 7 
and use this to relate the expressions for D\ X L and D 22 L = f V' 2 + V 2 . 

In particular, we begin with the equation (4.9) and then start by defining 
a particular primitive U of V. The procedure is in two steps, as we need to 
correctly choose the constant for the primitive. So first we set 

(4.11) u(y) = f ' V{s)ds 

J a 

Note that we need to check the well-definedness of u on 7, as it is a closed 
loop; on the other hand it is enough to check that the period u(2ti) — u(0) = 
J 7 V vanishes (here using the obvious notation for a pair of endpoints for 
the loop). 

For convenience in the sequel, set 

Im $ 

(4.12) F = —. 

9o 

so that equation (4.9) becomes 

(4.13) Vyy ~ V = ~Ty 

Then, for well-definedness of u, we observe that (letting subscripts indicate 
differentiation in the variable) 



n(2vr) - u(0) = I V 

J *y 

= f(V v + Hydy 



n 

(4.14) = 0. 

Thus u (and u + c, for any constant c) is well-defined along 7. 
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Next begin again with the equation 

(4.15) Vyy-V= ~T y: 

and then note that 

(Uyy -U + J r )y = Vyy-V+J r y 

(4.16) = 0. 

Thus we have that u yy — u + T = Co, where cq is a constant. In particular 
if we set 

(4.17) U = u + c , 
to be another primitive of V, then 

(4.18) Uyy-U = ~T. 

The point of all of this is that the positive part of the second variation of 
length integral (4.1) will turn out to be the energy of U, while the negative 
part will once again be the L 2 norm of T along £q (cancelled out by a term 
in the metric variation contribution D 2 ±L). 

We compute the contribution —D\ 2 L from the second variation of length 
along the surface through 



(4.19) 

-D 2 22 L = - [ V 2 + V 2 



10 

= / VyyV - V 2 by parts 

'70 



= / (Vyy ~ V)V 

= I {-Fy)V by (4.9) 

J 70 

= / TV y by parts 

•^70 

= / TU y y from the definition of U as a primitive of V 
= I T{U-F) from (4.18) 

•^70 

= - f T 2 + f XJT 

J 7(, J 7(, 

= - ( F 2 + ( U{-(U yy - U)} from (4.18) 

J 70 J 70 

(4.20) = - I F 2 + f U 2 + U 2 by parts. 

J -vn J Tn 
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Combining this last equation with (3.5) and (4.1), we find that 
d 2 

-j^L(g t ,j t ) = L>i 1 L{g ,'yo)\g,g\ - D^ 2 L(g ,-f )[i][i] 



T 2 



70 



2(A - 2) 



-i 



(4.21) 



2(A- 2)- 1 



|$| 2 



J lo 



Vy + V 2 



21 



(4.22) 



f T 2 

-- I T 2 + f U 2 + U 2 from (4.2) 

J 7 J 7 

f -2(A-2)-^+ f U 2 + U 2 . 

J 7o % J 70 



In summary, the Weil-Petersson Hessian of length can be expessed as the 
sum of two integrals along the curve, each of which has a positive function as 
an integrand. The first integrand is the restriction to the curve of a solution 
of a differential equation on the surface, and the second is the energy density 
of a solution of a differential equation along the curve. 

We record this formula as a theorem, extending Theorem 1.1 from the in- 
troduction. To set the notation, let [7] be the free homotopy class of a closed 
curve (simple or not) on the surface, and F(t) a Weil-Petersson geodesic arc; 
the class [7] is represented by the r(t)-geodesic 74. The tangent vector to 
Teichmiiller space at T(0) is given by a harmonic Beltrami differential, say 
— . Let — denote a second harmonic Beltrami differential on r(0). 

So So v ' 

Let J7* and denote the respective solutions to the ordinary differential 
equations (see (1.1)) 

Im$ 

(4.23) 
and 
(4.24) 



Uyy-U= - 



90 



U, 



yy 



U 



Im* 

5o 



This is enough terminology so that we may summarize our discussion as 



Theorem 4.1. Along the Weil-Petersson geodesic arc T(t), the second vari- 
ation ^1 of the length £(t) = L(T(t), [7]) is given by 
d 2 r 1*12 



(4.25) £2*®= I -2(A-2)- 1 l^ds+ / [U*] 2 + [U*] 2 ds. 

aZ J 70 9o J 70 

More generally, the Weil-Petersson Hessian HessL[^, J^] is given by 

(4.26) HessL[-,-}= [ -2(A - 2)- 1 ^^ds + f U?U* + U*U*ds 

90 90 J™ 9n An y y 
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Proof. The solutions 17* and 17* to (4.23) and (4.24) are unique, and the 
equations are linear in the unknown and the parameters <I> and Thus the 
unique solution to 

(4.27) Un . U = Jj^±ll 
satisfies 

(4.28) [/*+* = £/* + 17*. 

Then a straightforward polarization of (1.2), together with our understand- 
ing (4.28), yields (4.26). □ 

Remark. In terms of our previous notation for the Beltrami differential [i = 
J^, the equation (4.23) takes the form 

(4.29) U yy - U = Im fi. 

4.3. A geometric kernel representation. The second term of equation 
(1.2) is expressed as the energy of the solution of a differential equation. 
We wish to provide a more geometric interpretation; not only do we hope 
that this version is more appealing on its own, but it will be important in 
section 9 where we treat the Thurston metric via the Hessian of the length 
function. 

To begin, note that the second term of (1.2) (=(4.25)) may be written 
f U 2 y + U 2 ds = - f U{(U yy - U)ds 

J 7n J 7n 



'70 ^70 

U{s)F{s)ds 



I 



in the notation where 7 r (s) = Im = — Im /i, since 

U vv - U = -T 



by (4.18). It is well-known that we can represent the solution U(s) to (4.18) 
by 

U(s) = - j F(t)K(s,t)dt 
for kernels K(s, t) which satisfy 

(4.30) ^K(s,t)-K(s,t) = 6 s (t). 

It is easy to guess the solution to (4.30) using that the solution to ^ I K(s, t) — 
K(s,t) = are linear combinations of sinh(t) and cosh(i). (Indeed, if we 
represent 70 as the interval [—7/2,7/2] with endpoints identified, set so = 
±7/2 to be an endpoint, and look to solve (4.30) on that interval, then it is 
evident that setting K(s,t) = — cosh(t) is correct up to an easily computed 
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multiplicative constant.) In general, for 7 parametrized by an interval of 
length L (so that we may choose \t — s\ < L/2, we have that 



K(s,t) = 



1 cosh(s-t-L/2) , 

2 sinh(^) ' 1 ^ b 



1 cosh(t-s-L/2) 

2 s inh(-|) 



t > S 



solves (4.30) (where we require that \t — s\ < L/2). 

Of course, the variables s and t parametrize the curve 70 with respect to 
arclength, and so, for \s — t\ < L/2, we have \s — t\ = <i(7o(s), 7o(i))- Thus 
K(s, t) admits the description in terms of p = 7o(s), q = jo(t) as 

This leads to the representation 

d 2 f l^ 12 

u -£(t) = I -2(A-2)~ 1 i ^ 

70 9o J"/o 



dt 2 



(4.32) 
(4.33) 



l(t) = f -2{A-2)~ 1 ^- + f Ul(s) + U 2 (s)ds 

J 7o 9q J 70 

= I -2(A-2)- x ^-- f U{s)lm l i(s)ds 
J 70 <?o J 70 

= f -2(A - 2)~ 1 ^- - f Im/x(s) / K(s,t)lm l i(t)dtds 

J 70 <?0 J 70 ■/ 70 

= / -2(A-2)" 1 ^-- / / Im//(s)tf(s,t)Im/i(t)dtds 
V 70 <?o JJ 70 x 70 

= I -2(A-2)- 1 ^-- ff Im n{p)K{p,q)Imn{q)ds{p)ds{q) 
J 70 9o J J 70 x 70 

r 1 i$i 2 

= / -2(A-2)-^ 



+ _ . \tl\ ff lmfi(p)[cosh(d(p,q) - ^)]Im fi(q)ds(p)ds(q). 
in ( 9") J J70X70 ^ 



2 sinh( ^) j ^70x70 
where ds(p) and ds(q) refer to arclength measure 



Part 2. Extensions and Applications of the formula for the 
Hessian. 

5. A LOWER BOUND EXPRESSED IN TERMS OF POINTWISE QUANTITIES. 

We claim 

Lemma 5.1. Let v = 1/3^-. Then v is a sub solution of (A — 2) / = 
and in particular < v < — 2A — 2)~ 1 (-^p-). 
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We begin by noting that the curvature of a metric expressed as G|(te| 2 is 
given by 

K{G\dz\ 2 ) = -II A log G 

where A = d 2 + d 2 . 

Then using that K(go\dz\ 2 ) = — 1 and that |$o| \dz\ 2 is a fiat metric 
with concentrated (Dirac function type) curvature singularities at the ze- 
roes <I> _1 (0) of 3>, we see that 

|$|2 

A log^- = A log|$| 2 - A log^ 
% 

= -A\<S>\K(\<$>\\dz\ 2 )+A 9() K(g ) 

= 4|$| K5 p deg p <5> - 4g , 

pe*- 1 ^) 

where 5 P indicates a delta function at p. On the other hand, using that 
Ao log F = — ^ V ^P , we see we may write 

,$,2 A ^ |Vn^' 2 



Aolog^f = — 

9 <> llg" 



( ! 

Putting the last two of these equations together yields 




$|2 |V 



vr5 p deg p $ -4^ + 



m 2 ) i2 

~9o ) 



In particular writing A = for the go-Laplace Beltrami operator on 

S, and noting the vanishing of the first term on the right hand side, we 
conclude that 

|$|2 |$|2 
A l -4- > -4 1 



2 — 2 ' 

ffo % 

We are of course interested in the operator A — 2, so we note the obvious 
implication that 

|$|2 |$|2 
(A-2)^>-6^ 
% % 

so that v = l/S^r is a subsolution for the equation (A — 2)f = —2-^-. 

It is obvious that v = > 0, and if / satisfies (A — 2)/ = — 

then A(/ — v) < 2(f — v) and so the minimum principle guarantees that 
at a minimum of (f — v), we have / — v > 0; hence / — v > everywhere, 
concluding the proof of the lemma. □ 
Combining Lemma 5.1 with Theorem 1.1 we obtain 
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Corollary 5.2. Let <I> G QD(go) be a holomorphic quadratic differential in 
(S, (70) a^rf Zei denote a Weil-Petersson geodesic arc with initial tangent 
vector given by the harmonic Beltrami differential Qg^ 1 ■ Let £{t) denote the 
geodesic length of a representative 7t of a curve class [7] on S. Then, for 70 
the geodesic represented of [7] on r(0), we have 



We will apply this estimate in the section 8. 

6. The Second Variation of the length of an Arc 

In this section, we adapt our derivation to the case where S is a surface 
of finite genus with a finite number of punctures, and we are interested in 
the variation of length of an arc a that runs between two of the punctures 
(or a puncture itself). Naturally, the length of such an arc is infinite, so 
we will be discussing the variation of some regularization of its length; nev- 
ertheless, all of the basic considerations will extend to this case with only 
minor modifications. 

6.1. Notation and preliminaries. Let at be the geodesic on (S,gt) that 
connects punctures p and q in a fixed homotopy class (rel p and q) . Consider 
a sequence of points p n , q n G ao with p n — > p and q n — > q and let a.t, n denote 
the finite length ^-geodesic arc connecting p n to q n which is homotopic (rel 

Pn, Qn) tO 

cto,n C ao- Our plan is to derive a formula for 



show that the limit exists and is independent of the choice of the sequence 

{Pn,Qn}- 

We learned while preparing this manuscript that Wolpert [Wol07] recently 
treated the analogous case of finite length arcs between horocycles. 

It is easy to check that the formal preliminaries remain the same as in 
the derivation of (2.4), and so we conclude 

d 2 

(6-1) —L(g t ,a t , n ) = D^L^o, a , n )[g, g] - D% 2 L(g , a , n )[ao,n, a 0>n ]. 

6.2. The Second Variation of Arclength of ao, n in gt- As in the case of 
a closed curve, the first term in (6.1) is relatively straightforward to compute; 
the only new issue to consider is the dependence of the term on the choice 
of endpoints p n , q n of a ti „. Indeed, exactly as in the derivation of (3.5) , we 
formally compute 




□ 



—L(g t ,a t:n ), 



(6.2) 




where the principal issue is to determine the meaning of (A — 2) 
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6.3. Variations of metrics of finite area. The basic point here in under- 
standing — 2(A — 2)~ 1 i^- is to construe it as H for the family of pullback 

metrics gt in (3.1). As these maps id : (S,go) — > (S,gt) are harmonic, we 
can apply some results from the theory of harmonic maps between cusped 
hyperbolic surfaces. 

In this direction, results in [Wol91] (Theorem 5.1) and of Lohkamp (see 
the remark after Theorem 4 in [Loh91], especially with Lemma 12 informed 
by Proposition 3.13 in [Wol91]) proved that H(t) G C k ' a (S,g ) was analytic 
in t on the compactified surface M. In particular H is bounded. 

Indeed, we can easily show from this that TL = 0( ^u a ) = 0(y~ a ) 

for some a G (0, 1) as r — > 0. The basic elements of this argument is 
that ^Tiys = U a f° r some a G (0, 1) is a supersolution of the equation 

(A — 2)7Y = — 2||<I>|| 2 on the cusp, as well as the point that the kernel of 
(A — 2) on a half-infinite cylinder C = {Imz > 1, | Rez| < 1/2} (with the 
standard identifications) is spanned by the pair of functions ki(z) = y 2 
and k>2{z) = y^ 1 . With that background, consider on the finite cylinder 
{1 < Imz < y n }, a function Hj (z) of the form Hj (z) = Coy~ a +C\y~~ 1 +ejy 2 . 

Then for appropriate choices of ej — ► 0, we find that Hj(z) majorizes 7i; 
letting j — > oo and tj — > while Cq and C\ stay bounded (as the boundary 
values H(z) for {Im z = 1} are fixed independently of j) allows us to conclude 
that H(z) decays like C y- a + Ciy' 1 . Thus H(z) = 0(y- a ) = 0( ^ l)a ). 

Looking ahead to the final form of the second variation of length, it is 
worth recording the 

Proposition 6.1. Let p n — > p and q n — > q, and let ao,n be the geodesic arc 
connecting p n to q n as in the introduction to the section. Then 

[ -2{A - 2)- 1 ^-ds = I hfLds 
converges as n — > oo. 

Proof. In the upper half plane coordinates, we have that for each end of an 

(6.3) f Hds= f yn H^L= 0{-^)— = 0{l) as n ^ oo. 
Ja n J a y Ja y" y 

The proposition then follows from the integrals being positive. □ 

6.4. The Jacobi Field for an Arc. The next term we must address is the 
second term in (6.1). The variational vector field, defined geometrically, also 
satisfies (4.9). Here, of course the variational field V n depends on n, as it is 
defined in terms of at jn ; our notation is meant to reflect that. The form of 
the second variation of arclength is also unchanged at f V 2 , + V 2 , as the 
boundary term vanishes once we require the family of curves a,t, n to have p n 
and q n as endpoints, independently of t. 
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The primitive U n = f% V n {s)ds exists as before — in fact, in the setting 
of an open arc, there is not a well-defmedness issue to check, though we still 
have to adjust by a constant. In particular, note that for U n = JJ V n {s)ds 

(so that U n = P V^'(s) - (Im fi)'(s)ds), we have U% -U„ = Im/i + (V^(p„) - 
Im//(p„)). So set J7„ = V n (s)ds + (V^(pn) - Im/x(p n )). Then U n solves 
the (boundary value) problem 

(6.4) [/; , -C/ n = Im/i 

^n(Pn) = V n {Pn) ~ Im/i(p n ) 

^n(gn) = Vl(q n ) - Imn(q n ). 

The last boundary condition follows after applying the fundamental theorem 
of calculus to 



rqn 

U n {q n ) = j V n (s)ds + V^ipn) - Im/x(p n ) 

J Vn 



Pn 

V^{s) - (1m Li)'(s)ds + VM - Im fi(p n ) 

Pn 

We will soon show that these boundary values, though non-zero, tend to 
zero and have no effect on the limiting relation. 

At this stage, it is useful to write a kernel representation for V n . Let 
L n = d(p n ,q n ), and parametrize the arc «o,n by [— %S^f]- It will not 
ultimately affect the results if in these coordinates the midpoints of [p„,<?„] 
(represented in our coordinates by the origin) do not remain in a compact 
set. 

The kernel for the operator J^- — 1 on the segment [— , 4f] is given by 



sinh(^+s) sinh(^--j/) 



(6.5) K n (y,s) 



sinh(L„) 

sinh(^-£)sinh(^+2/) L„ / „, / „ ^ L„ 



< S < y < 



< V < S < -Hf-. 



sinh(L„) ' 2 — » — ° — 2 

This gives the representations 

(6.6) V n (y) = r K n (y,s)(Im f i)'(s)ds 

JPn 

and 

(6.7) U n (y) = K n (y,s)(lmfi)(s)ds + a„cosh(y) +6„sinh(y), 

where a n and b n are chosen to satisfy the boundary conditions in (6.4). 

We will estimate the asymptotics of these boundary terms in the next 
subsection, but we display the preliminaries for that analysis here, in the 
present context of integral formulas for the relevant geometric objects. 
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One of the terms in the boundary condition is given by at the corre- 
sponding boundary point. For example, V^(p n ) is given by 



(6.8) *<*.>-£"| 
We then compute that 



K n (y,s)(Imn)'(s)ds. 

V=Pn 



(6.9) j-K n (y,s)=- 



sinh{Ln) , 2 - s - y - — 

sinh(^p- a )cosh(^p.+y) ^ ^ 

sinh(L n ) ' 2 - « - 6 - 2 



Combining (6.8) and (6.9) yields the representation 

sinhf-^ - s) 

(6.10) <{p n ) = - I .\ 2 (r , ' (Im »)'(a n (s))ds. 

J_Lj L smh(L n ) 

The final matter is the analogue of the derivation of (4.2) from the varia- 
tion of arclength. This involves three integrations by parts, and so we need 
to consider the boundary terms from each integration; however, all of the 
terms have either V n or V n ,y = V n as a factor, and so all vanish at the 
endpoints p n and q n . 

We conclude that 

(6.11) -Dl 2 L(g , a 0) n)(do,n, <*o,n) = - / + [ U^ s + U 2 



<*0,n J a 0,i 

in the notation of (4.12). 

We combine formulae (6.2) and (6.11) for the derivatives DnL(go, ao, n ) 
and £>22-^(5o,Q;o,n) to find 

(6.12) ^ L (g t ,a t , n ) = I -2(A - 2)~^ + / f/ n 2 s + E/>. 

6.5. Passage to the limit. Naturally, we are interested in taking the limit 
of (6.12) as n — ► oo. That the first term converges in the content of Propo- 
sition 6.1. 

For the second term, we need to estimate the asymptotics of U n ,y = Vn 
and U n . We claim 

Proposition 6.2. The fields V n converge to a field V defined on the entire 
arc «o- The primitives U n may be chosen so that, not only does U n converge 
to a field U with U s = V , but also 



(6.13) f U 2 >s + U 2 ^ [ 



U 2 + U 2 < oo. 



Here U s and U n>s indicate derivatives of U and U n with respect to the ar- 
clength parameter s. 

Proof. It is useful to begin with an observation. 
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Lemma 6.3. Near an end of (S,go) with coordinates from {\z\ < 1}, we 
have M = 0(r(logi) 2 ) and = 0(r(log±) 3 ). 

Proof. In the coordinate disk, the holomorphic quadratic differential = 
| + h.o.t, while ^ = r 2 (logi) 2 and ^ = rlog^Jp along radial geodesies. 
The estimates are then immediate. □ 

Recall the terms a n and b n in (6.7) that adjust the solution U n for the 
inhomogeneous boundary conditions. Our next goal is to prove that these 
adjustments are asymptotically inconsequential. 

Lemma 6.4. a n ,b n = o(e~ c ~^). 

Proof. The goal is to show that the boundary conditions given in (6.4) are 
small. In that case, since we would have a n cosh(y) + b n sinh(y) = o(l) for 
y = ±-| from (6.7), we easily see the statement of the lemma. 

There are two types of terms in the expression (6.4): those given by 
|(Im/i)(p„)| and \(lm fi)(q n )\, and those given by V^(p n ) and V^(g n ). We 
will treat them separately. 

We first claim that \(Im/j,)(p n )\ and \(lm /j,)(q n )\ decay as n — > oo. This 
of course easily follows from Lemma 6.3. 

It is only slightly more difficult to use (6.10) to show that V^(p n ) and 
V' n {q n ) decay to zero as in oo. To see this, note that for any fixed y* , we 

have J y * ^si^L )^ ds = 0(eT~%~~ y ), and that |(Im/z)'| is bounded on all of 
the geodesic arc ceo, while vanishing into the cusp by Lemma 6.3. (Note also 
that J \ n sinh(L ) ds < 1.) These preliminaries are enough to estimate 

Vn(p n ) (with an analogous argument for V^(q n )). 
We begin from (6.10) with 

/y* sinhf^ — si 
l Jjf 1 > \(Im»y(a n (s))\ds 

< max lm/i / — — . — ds 

iPn,V] J-*p smh(L n ) 

"o J y , smh(L n ) 

Then, for e small, use Lemma 6.3 to pick y* so that max^^.j |(Im//)'| < e. 
Then the first term is bounded by e while the second term is bounded 
by maXoiQ | (Im//) / |0(e~~2 L_J/ ). Letting L n — > oo and — f- < y* — ► — oo 
somewhat more slowly than —k^- — > — oo (e.g. y* = —ty) shows that 

K(Pn)^0. ' " □ 



20 



M. WOLF 



With these preliminaries, we easily conclude the proof of the proposition. 
As L n — ► oo, we find that for fixed y, the kernels K n (y, s) limit on 



K(y,s) 



s <y 



2 ' 

p y — s 

y < s 



which we can write succinctly as 

As both (Im / u)(s) and (Im/x)'(s) decay while K n (y,s) converge to an in- 
tegrable function, we see that the formulas (6.6) and (6.7) show uniform 
convergence of V n to a well-defined variation field V, as well as primitives 
U n to a well-defined (primitive) function U. As U' n = V n , the uniform con- 
vergence of {[/„} and {V n } show that U' = V. Using Lemma 6.4, we obtain 
the formula 

U(y) = f K(y, S )(lm^(s)ds. 
(6.14) = / 9 le- d(y ' s) {lm^(s)ds 

Next, it is easy to see that U{y) — ► as either y — ► p or y — > q. Mimicking 
the argument that showed that V^(p n ) — > 0, choose y* close enough to p 
so that max[ pj/ *] < |. Then choose y even deeper into the cusp so that 

%,y*)>-log(;^^)Then 

|£%))<| f V K(y,s)(lm^( S )ds\ + \ P K(y, s)(lm ^)( S )ds\ 

Jp Jy* 

oo 



/-i/ /»oo 1 

<-/ |if(y,s)|ds + max|/x| / -e-^ds 
2 J p Q o Vy* 2 



< J + Je-^^maxI//! 

2 2 ao 

< e 

as desired. 

Finally, we address the finiteness of the energy of U. Of course, the metric 
go near the cusp point p (or g) may be expressed as go = |z| _2 (log |z|)~ 2 |(fc| 2 , 
and so the estimates for \x may be written as \n(y)\ = 0{e y ~ e ). It is then 
an easy estimate that 

I U' 2 + U 2 = ]- 1 1 e~ ls - yl Im /i(s) Im ^(y)dsdy < oo 

and is the limit of f U' 2 + U 2 . 

This concludes the proof of the Proposition. □ 
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We summarize our discussion in this section with formulae for the second 
variations of an open arc a analogous to those in Theorem 1.1 for the second 
variations of length of a simple closed curve 7. 

Theorem 6.5. Along the Weil-Petersson geodesic arc T(t), the second vari- 
ation jjpi of the r (t) -length £(t) = L(T(t),a) of a (class of an) arc a is given 
by the (convergent) expression 

(6.15) 

4*(t) = / "(A - 2)-^ds + / [tf] 2 + [U^ds 

al Ja 9q Ja 

(6.16) = [ -(A-2)- 1 ^^ds + - [[ e-^ s ' y hmfi(s)lmfi(y)dsdy. 

J a 9q ^ J J aQ xao 

More generally, the Weil-Petersson Hessian HessL[^,^£3Lj i s gi ven by 
the (convergent) expression 



(6.17) HessL[-,-}= [ -(A - 2)' 1 2Bb ** ds + I U*U* + U*U*ds 

90 go Ja 9o Ja " " 

f / a 1 2 Re < I ) \& 

(6.18) = / -(A-2)" 1 2 ds 

Ja 9o 

+ - [I e _(i ( s 'f) Im/i(s) Imu(y)dsdy, 

2 J Janxan 



where /J, = ^ and v = are the harmonic Beltrami differential representa- 
tives of two tangent directions at [go] . 

Here satisfies the equation (1.1) along the arc 7, a condition which 
forces ?7* — ► along 7 as it tends to the punctures p and q. 

7. Convexity for Laminations 

We have already seen in Theorem 1.1 and Corollary 5.2 that if 7 is a 
simple closed curve, then on a Weil-Petersson ray T = T(t), we have 

(7.1) ^ i{T {t))> l -fm\ 2 d S 

where $(t) is the holomorphic quadratic differential tangent to V at T(t), 
and || $|| = In this section, we extend that result to prove the 

Proposition 7.1. Let V = T(t) be a Weil-Petersson ray and A a measured 
lamination on S. Then 

(7.2) £<»<rM)>ij(W* 

where Q t is the holomorphic quadratic differential tangent to T at T(t). 
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Our definition of J A ||<fr|| 2 ds is straightforward and parallels the definition 
of length of the lamination A. See [BonOl] for a background discussion. In 
particular, a measured lamination A is defined as a measure X(k) = J k dX 
on transverse arcs k. Choose arcs k±, . . . ,kj which are transverse to A and 
construct flow boxes {Fi} for A bounded by the kj and parallel to A. 

Then if G\ is the geodesic lamination underlying A, then A — Ujkj is a 
(possibly infinite) collection of finite length arcs. The length of A is then the 
integral, with respect to the transverse measure dX of A, of the lengths of 
the finite arc components. More precisely, we lift each of these components 
A to A a C T 1 M, endow A a with the natural arclength measure ds, and then 
integrate the product to get 

£{X) =: [ ds=: ft dsdX(a). 

In order to define f x \\<&\\ 2 ds, we proceed analogously, except that we note 
that the function H^H 2 on S then naturally defines a measure H^ptis on 
T l M. In other words, we set 



f \\<S>\\ 2 ds = [[ \m 2 dsdX(a) 



Proof of Proposition 7.1. Let 7 n be a sequence of simple closed curves con- 
verging to A. The idea is to apply (7.1) to 7„ and then take a limit in n to 
find (7.2). 

Now £ 7 and £\ are real analytic functions on T [Ker85], and thus since 
£ ln — > £\, so does Jp-^ 7n — ► j^i^x- Thus the left-hand sides of (7.1) converge 
to the left-hand side of (7.2). 

For the right-hand side, the argument is virtually tautological. We first 
note that the arclength measure ||$|| ds is continuous on T^M. Then, choose 
n sufficiently large so that the flow boxes {Fi} described above also serve 
as flow boxes for j n . Then the definition of convergence in M.C then easily 
implies that the right-hand side of (7.1) converges to the right-hand side of 
(7.2). □ 

8. Convexity of £^ and upper bounds on the Hessian 
8.1. The function lh. 

8.1.1. The function £3. In light of (7.1), we can quickly refine the basic 

convexity result for the length £ to a convexity result for a concave function 

2 

of £, namely £3. We will see in the next subsection that this is not sharp, 
but at this stage it is elementary. 

To begin, recall from (3.7) that the first variation of length may be ex- 
pressed as 



d f Re<I> 
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We can then estimate this derivative as 



\±UT(t))\ < I \ — \ds 
at J l0 go 



70 

cf>l 2 , 1 I 

70 



<(/ t^dsynl 

J 70 % 



Squaring and combining with (7.1) yields 

(8.i) 4U( r (*))> i(4^(r(*))) : 



V „ _ 

dE 1 



We compute J^3(r t ) and substitute in the above inequality to conclude 
the 



2 

Corollary 8.1. The function £/* is Weil-Petersson convex on Teichmiiller 
space. 

8.1.2. The length of an annulus. In this subsection, we compute in two ways 
the second variation of the length of the core geodesic in a hyperbolic annu- 
lus; the basic result is well-known (see Example 3.6 in [Wol08]). 

let C{£) denote a complete hyperbolic annulus whose core geodesic has 
length I. This cylinder may be parametrized as 

CTO = [-^]x[o,i], 

where top and bottom edges are identified, and we consider the metric cylin- 
der as equipped with the hyperbolic metric g = dsf = £ 2 esc 2 £x\dz\ 2 . 

Consider rotationally the harmonic map R(t) : C(£) — > C(£ + t) which 
does not twist the boundary; in other words, this map may be expressed in 
coordinates as R(t) = u{t) + iv{t) where u{t)(z) = u{t)(x) and v(t)(z) = 
y. Now the rotationally invariant holomorphic quadratic differentials on a 
cylinder have a particularly simple form: we may write one as <1> = cdz 2 
in the complex coordinates above. Thus, taking one of these as a Hopf 
differential for our map, and using that the holomorphic energy TC, the 
Beltrami differential v and the Hopf differential <!> may be related by = 
g(£o)Tti>, we conclude that 

, 1 + v I 2 + c cos 2 £x 
\ — v £ 2 — c cos 2 £x ' 
If c(t) is the factor so that the Hopf differential is parametrizing a family of 
harmonic maps R(t) : C{£) — > C(£ + t) whose targets are progressing through 
Teichmiiller space T(C) at unit Weil-Petersson speed, then two conditions 
hold: (i) the choice of c(t) provides for c(t)dz 2 to the Hopf differential for 
the map R(t), and (ii) || j|$(£)|| wp = 1- 

Now, for w(t) to have image C(£ + t), we must have the boundary of C{£) 
map to the boundary of C(£ + t), i.e. 

TT [% 1 + V l 2 + CCOS 2 £ X 

7ZT7, r = u — = / u (x )dx + u U = / = Tnrdx. 

2(£ + t) y 2£> J y ' v ' J l-i/ £ 2 -ccos 2 £x 
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Upon differentiating in in t and finding the resulting elementary integrals, 
we obtain c = — I. Thus 

\\-QgWwp = Wcg^Wwp = \\(-£)£~ 2 cos 2 £x\\ W p = — 

after another explicit integration. Thus || J^||jyp = "t - and so dsyyp = ^d£ 2 

on the Teichmiiller space T(C). This implies that on this space, ^ = tt~ 1 £^ 

and so I = (2n)~ 2 s 2 . Thus the length £ of the core geodesic satisfies that 
is convex, but not convex to any lower power. 

Remark. Alternatively, we may use the formulas (1.2) and (3.7) to ana- 
lytically find the same result. In that case, if we set = cdz 2 , then 
i = f — 2(A - 2)~ 1 {c 2 g~ 4 }. However, the equation Au - 2u = c 2 g~ 4 
reduces in this case to an ordinary differential equation, whose solution we 
can require to be bounded on the boundary. It is then elementary (see the 
analogous analysis in (8.18)ff using the method of variations of parameters) 
to find an exact expression for m on 70. We then compare with the expression 
in (3.7) for the first derivative to to obtain the half-power convexity result 
above: this confirms, at least in this very simple case, the formula (1.2). 

8.1.3. £2 is convex. We now offer a geometric argument of Wolpert's recent 
1 

result [W0IO8] that is Weil-Petersson convex. 

Comparison of lifted harmonic map to rotationally invariant map. The es- 
sential point is best understood in the setting of the annular covers (C,g~t) 
of the family of surfaces (S,gt). Consider the harmonic maps wt : (S,go) — > 
(S,gt) and their lifts wt '■ (S,go) — > (S,gt). These lifts are in the homotopy 
class of the rotationally invariant harmonic map R t : C(£q) — ► C(£t), where 
ty(gt) = £t- Now wt is conformal only at the (isolated) zeroes of the Hopf dif- 
ferential, and so, off of small neighborhoods of the zeroes of the lifted Hopf 
differential, the harmonic map Wt has quasi-isometric constant uniformly 
bounded away from 1. By contrast, one can either compute or reason geo- 
metrically that the rotationally invariant harmonic map Rt : C(£q) — > C(£t) 
has quasi-isometric constant tending uniformly to 1 as one leaves compacta 
in C(£q): the image curves are growing exponentially in length, so energy 
efficiency requires the map to be increasingly close to an isometry as one 
leaves compact sets. 

Let 7i R (t) be the holomorphic energy (see equations (3.1)-(3.3)) of R(t) 
and H(t) be the holomorphic energy of Wf Since both R(t) and wt are the 
identity when t = 0, and since, as we have just seen, R(t) is asymptotically 
an isometry while wt is boundedly away from being the identity off small 
sets, then we find 

H(t) > H R (t) 

outside some large compact set (at least away from small neighborhoods 
of the zeroes of the lift $ of <E>). In particular, parametrizing C(£) as in 
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subsection 8.1.2, we see that 

(8.2) f H(t) > [ H R (t). 

A comparison of ODEs. The rest of the proof follows by applying other 
inequalities that reflectthat R(t) is a harmonic map of lower (regularized 
in some way) energy that w t . In particular, consider the Fourier expansion 
$ = J2^n(x)e 2nmy of the quadratic differential <E> (where the map wt has 
Hopf differential $). 
Because 

f Re<l? 

£= / ^dy = £Reb 

Jx=o 9 

we know that the Hopf differential <E> R for the rotationally invariant map 
R(t) must be & R = Hebodz 2 on C{£): this is because the targets C(£ + t) 
agree for the two maps wt and R(t), and hence the change in core-curve 
length is the same. 

The upshot is that, for an arbitrary constant curvature circle {x = £}, we 
have 

(8.3) f \^ r ds = g-i{x)Y J K\ 2 >g-'i{x){Reb ) 2 = I ^-ds. 
Jx=£ 9 ' ' Jx=( 9 

Of course, we know from formula (1.2) that 

f i l$l 2 1 f •• 

£> -2(A-2) _1 ^-ds = - / Hds 
J 7 9 2 7 7Q 

To estimate this last integral, let u be the solution of 

-2|$| 2 

A„u - 2u= ^ — . 

g 2 

If we were to integrate this equation along the vertical parameter curves 
{x = const}, we would obtain an ordinary differential equation for the 
function f udy = f ^Hdy in the single variable x G (— ^, ^), i.e. 




Of course, a similar equation holds for the integrals f hH R dy and f '® J dy 
associated to the rotationally invariant map. Indeed, inequality (8.2) says 
that there are boundary points x = ±^ =F 5 at which J x ^H R dy < J x \Hdy; 
moreover, inequality (8.3) asserts that on the interval (—^e, ^s), the right- 
hand-side of (8.1.3) is less in the lifted case than in the rotationally invariant 
case. 
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The upshot is that the comparison principle for ordinary differential equa- 
tions implies that 

|2 



f _2(A-2)- 1 %dy= f\H R dy 
Jx 9 J X z 



>f X 



In particular, specializing to the curve {x = 0}, and recalling the implication 
above of (1.2), we find 



j-f y Jj y <- 



Here the last inequality is inherited from the rotationally invariant case: 
recall that our choice that Re f $> = f 3> R implies that the infinitesimal 

J 7o J 7o r 

change of lengths agree between the lifted and rotationally invariant maps. 
We conclude 

Corollary 8.2. (Wolpert [WolOS]^ The function 1^ is Weil-Petersson con- 
vex in the Teichmiiller space T(S). 

8.2. A general upper bound for the Hessian. We have already seen a 
lower bound for the Hessian of length in Corollary 5.2 and Proposition 7.1. 
In this passage, we note an easy upper bound as well. 
We begin by noting that if 

|$|2 

(A - 2)h = -2 l -4r- 
9o 

on a surface 5, then the maximum principle implies 

|$|2 

(8-4) /*< H^Hoo 

9o 

where the right hand side is the maximum of the function J — 3- on S. In the 

formula (1.2), this will estimate the first term. 

To estimate the second term, we consider equation (4.18) (combined with 
(4.12)) 

Im<I> 

(8.5) u U = - 

9o 

The maximum principle then implies that 

,Im$ , 



(8.6) U < max 



7 9o 
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Thus the second term in formula (1.2) is estimated as 



8.3. Estimates for the Weil-Petersson connection near the com- 
pactification divisor. In this passage we refine the method above to es- 
timate the Weil-Petersson connection on a codimension two distribution 
V C TM of the tangent bundle near the Deligne-Mumford compactification 
divisor (i.e. for surfaces with small injectivity radius) which is in some sense 
"parallel" to the tangent bundle of the compactification divisor. Roughly, we 
prove that V is quite flat, in the sense that for X, Y G V, we will have that 
the normal component (V^^) 1 of (VxY) satisfies (Vx^) 1 = 0(£ 2 ). (Here 
£ signifies the length of the curve which vanishes on the nearby component 
of the compactification divisor. 

To state this precisely, we choose a simple closed curve 7 C S and a small 
number £ > 0; we consider the level set L^{£) of hyperbolic surfaces for which 
L(-,7) = £. The set L y (£) is a submanifold of the Teichmiiller space T of 
real codimension one, and it is orthogonal to the vector grad£ 7 , the Weil- 
Petersson gradient of £ 7 . Let J be the almost complex structure of T and 
consider the projection r G TL^(£) of Jgrad£ 7 into TL^it). Let V C TL 7 (£) 
denote the distribution of (dimT — 2)-planes in TL^(£) orthogonal to the 
span of grad^ 7 and r; note that V is Weil-Petersson orthogonal to both 
grad£ 7 and Jgrad^ 7 . 

Remark. One does not expect V to be integrable. In particular, one does 
not expect that Jgrad^ 7 is parallel to the Fenchel-Nielsen twist vector field. 
(See [Wol82].) Nevertheless, it will be a consequence of equation (8.10) of 
the early part of the next proof that, as £j — > 0, the distribution V converges 
to the tangent bundle TC 7 of the compactification divisor C 7 = {^ 7 = 0} of 
the augmented Teichmiiller space T. (Compare [W6191].) 

Now, for X G V and Y a section of V — > T, (i.e. a vector field on 
Lj(£)), we can consider the Weil-Petersson covariant derivative VxY. Of 
course, the vector VxY has components both in V and in the orthogonal 




after a substitution 



= L(max I 1) 2 . 

7 So 

We conclude, taking into account Corollary 5.2, that 
Corollary 8.3. 
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complement V of V; we focus here on the component of this vector in the 
orthogonal complement. 

Remark. It might be difficult to formulate general results on the full vector 
VjF in useful and incisive way. For example, if we were to "lift" a curve a 
from the compactification divisor C 7 to an almost parallel curve a tangent 
to V, then since Vo-d C TC 7 can be arbitrary, so might we expect V^d to 
have an arbitrary non-orthogonal component. 

Our main result in this section is 

Theorem 8.4. In the notation above, for X £ V and Y C V of unit norm, 
we have 

(8.9) (Vx*0 X =< VjF, V ± >= 0(4). 

In particular, < Vx^gradty >= 0{£^) and < VxY, Jgrad£ 7 >= 0(£^). 

Remark. Similar results were obtained recently by Wolpert [W6108], using 
his estimates on the Hessian. From the formula for the Hessian presented in 
(1.2), we see in the present derivation the elementary nature of the expan- 
sion of VjfF in £j. The flatness of order 0(£^) occurs because the (explicit) 
rotationally invariant even solution u(x) = xtan£x + l/£ of the Jacobi equa- 
tion has the following property: at the core geodesic of the cylinder, this 
function u(x) is smaller by a factor comparable to 0(£ 7 ) than it is on the 
boundary of the cylinder. Since the geodesic has length 0(£ 7 ), the dominant 
term of the integral of u(x) over the geodesic decays like 0(£^). 

Proof. We begin the proof of Theorem 8.4 with a preliminary proposition 
characterizing the quadratic differentials which represent elements of V. 

Proposition 8.5. Let X G T^ M ^T with X G V and <I> = &x be a holomor- 
phic quadratic differential on the surface (S,g) for which <&/g is a harmonic 
Beltrami differential representing X. Then for 7 the geodesic on S used to 
define L 7 (£), we have 

(8.10) / -ds = 0. 

J y 9 

Proof. Since X G V, we have < X, grad£ 7 >= 0. But from (3.7), we find 
that 

=< X,gradty > 
= X(Q 

f Re$ , 
= / ds. 

J 7 9 

Of course the other defining condition of X G V is that < X, J grad £^ >= 0. 
But since the Weil-Petersson metric is Kahler, this implies 

=< X, J grad £^ > 

= — < JX, grad £~f > . 
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But JX is represented by i<&/g, and so we derive as above that 

0= / ds, 

Jy 9 

proving the result. 

□ 

Remark. The covectors in T*C 7 are represented by holomorphic quadratic 
differentials which have simple poles at the nodes. On the other hand, a 
generic covector in T^T at an element S € C 7 has a second order pole with 
a non-vanishing residue at the node obtained by pinching 7. 

Continuation of the Proof of Theorem 8.4 The heart of the matter is a 
computation of Kess£ 1 (X, X), in particular to prove 

Lemma 8.6. For X C P as above, we have 

(8.11) HessZ 1 {X,X) = 0{l 2 1 ). 

To see that this is enough, note first that polarization will imply then 
that Hesst^X, Y) = O(^) for X,Y CP oi unit norm. But then 

- < V x y,grad£ 7 > = -(V X Y)£ 1 

= (XY — VxF)f 7 since X,Y are tangent to L 7 (£) 
= Hess£ 7 (X, Y) by definition 

(8.12) = 0(4). by Lemma 8.6 
Moreover 

< V X Y, Jgradty > = - < JV X Y, gr&dl^ > 

as J is a Weil-Petersson isometry 
= - < VxJY, grad£y > 

as Weil-Petersson is a Kahler metric 
= - < Vx-Z^gradty > 

for some Z G P as P is J- invariant, being orthogonal to a J- invariant 
subspacc of a Kahler manifold. Then < VxY, Jgrad£ 7 >= 0(£ 7 ) follows 
in the manner of (8.12). This concludes the proof of Theorem 8.4, pending 
the proof of the main lemma. □ 

Proof of Lemma 8. 6. The basic idea of the proof is to estimate the terms 
in the formula (1.2) for the Hessian, where we take X to be represented 
by a harmonic Beltrami differential [i = &/g, and we apply the features of 
X C P C TL 7 (£) to prove that those terms are small. In particular, because 
the length L(g, 7) is small, the geodesic 7 is embedded in a wide, thin collar. 
Then, by Proposition 8.5, we learn that |$| must decay rapidly towards the 
center of the collar. Those facts together are enough to conclude that each 
of the pair of terms in (1.2) is small. 
We carry out the plan in steps. 
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Step 0. The collar. Consider the collar C = [— \ sec -1 j, \ sec" 1 \] x [0, 1] 
with horizontal edges [— \ sec -1 j, \ sec -1 \] x {0,1} identified, equipped 
with the hyperbolic metric go = I 2 sec 2 ix\dz 2 \. This collar C embeds in 
a neighborhood of the geodesic 7, with {0} x [0, 1] mapping onto 7. 
Step 1. Decay of |$|. On this annular collar C, we may regard the quadratic 
differential $ as a function (or more formally, we divide <E> by the nonva- 
nishing holomorphic quadratic differential dz 2 to obtain a function in the 
quotient). Then, by the rotational invariance of the collar (or by working in 
Fermi coordinates), we see that g may be taken as constant along 7, and so 
the conditions in Proposition 8.5 imply that <1? has no period in the collar. 
Finally, we estimate boundary conditions. 

Since, for t small, every X <ZV C. TL 7 (£) can be approximated on com- 
pacta away from 7 by an integrable meromorphic quadratic differential on a 
noded Riemann surface in the compactification divisor C 7 , and the collection 
of such quadratic differentials of unit Weil-Petersson norm is compact, we 
see that we may take |<&| as bounded on the horocycle of length one on dC. 

Because <3? is holomorphic and hence harmonic, the vanishing of the period 
together with the fixed boundary conditions is enough to show that |$| 
decays rapidly on the interior of the cylinder. 

We can obtain an estimate of this decay through a Fourier analysis of 
3>. On the cylinder C, set $ = a n (x)e 2mny . Then the harmonicity of $ 
implies 



Now as ao(0) = J x=0 $ = by (8.10), we find that J^, <3? = along any cycle 
7* in C 7 homologous to 7, hence ao(x) = 0. 

Of course, j dC |$| 2 < Cq by our argument on limits above, and so 



= A$ = 5^(a£(s) - 47rn 2 a„(x))e : 



2-niny 



and in particular the equations 



o"(x) — 47rn 2 a n (x) = 0. 



(8.13) 




We note that 



(8.14) 



(Y,a 2 n (x))" = Y,Ka n + 2K) 2 



n^0 n^0 



= ]T87rVa 2 +2(a'J 2 by (8.3) 



n^0 




n^0 
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Thus by the maximum principle applied to the differential inequality (8.14) 
with boundary conditions (8.13), we have 

(8.15) / |*| 2 < Co^hVfrrso ;= ^ cosh ^ 

Jx=x cosh V8tv(j sec 1 \) 

Finally consider $(zo), for zq = xo + iyo, a fixed point of the collar. Note 
that on our parametrization of the collar, the balls of radius ^ inject into 
the parameter domain. As $ is harmonic, 



|*(*b)| < -I / $1 

* J By (zn) 



JBi(zo) 

" J\Ba(z-zo)<$ 
Rc zo + l 



= -=( / / |$| 2 «2 

V 71 " JRc2 -| A=t 

^ /.Rc2 + 5 j 

<-^(/ Dcoshv^vrtdt)^ by (8.15) 

V 71 " Jr C20 -| 

= Di (cosh \/8vrxo) 5 



We conclude that 

|$(z )| 2 < D 2 cosh-\/87rxo, 

where (using (8.15)) 

(8.16) D 2 = 0(e~^) 

in for £ small, justifying our remark about the decay of \&(z)\ into the 
collar. 

Step 2: the function U. There are two terms in the formula (1.2) for the 
Hessian; here we estimate the one involving the energy of the function U. 
In particular, we already know from (8.7) that 



f u' 2 + U 2 < £(max^f 
J 7 9 



< ig~ 2 \ x =oD 2 

= o(r 3 e^) 

using g\ x= o = £ 2 and (8.16). This term is then consistent with the statement 
of the lemma that Ress£(X, X) < 0(£ 2 ). 
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Step 3: the function (A — 2) We are left to estimate the second term 

in the expression (1.2) for Ress£(X, X), namely 



i l$l 2 
2(A -2)- 1 - 1 



/ ^ 9 2 



In particular, we need to estimate the solution u$ to the equation 

|$|2 



(8.17) (A - 2)u = , . 

evaluated on the core geodesic. 

We again estimate the solution to this by estimating it on the cylinder 



2 



C 7 , on which we have good control on ^j-. Of course we already know from 
the maximum principle (see (8.4)) that 

|$|2 

Uq < SUp ■ 



9 2 



Now the latter is bounded C 7 , and the complement of C 7 has bounded 
geometry on which f Mr ^ c ^-dA g < 1 (because $ is of unit Weil-Petersson 
norm). We conclude that there is a Co for which 

|$|2 



9 



2 lac. 



7 



<C 



Thus by the maximum principle, it is enough to estimate the solution u of 
the boundary value problem 

1 „ D 2 cosh V8t:x 
-u (x) — 2u(x) = = 

g g 2 

n(±-sec- 1 -) = C* 

Using that g = I 2 sec 2 £x, we rewrite the equation above as 

(8.18) u"{x) - It 2 sec 2 lxu(x) = D 3 cosh( v / 8vrx) cos 2 ix 

1 1 
^i-sec' 1 -) = C , 

where D3 = 0(£~ 2 e~ v ^ 7T / e ). The homogeneous equation 

u"(x) - 2£ 2 sec 2 lxu{x) = 

has the two solutions 

(8.19) u 1 {x)=taxi£x 

u 2 (x) = x tan £x + — . 
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Using these, one can solve for a particular solution of the form uq = u±vi 
U2V2 by elementary integrations, namely 



vi = J U2-D3 cosh(\/8vrx) cos 2 £x 
V2 = — J U1D3 cosh(v / 87rx) cos 2 £x. 



An asymptotic expansion shows that u\V\ + u 2 f 2 |,9c 7 = 0(1); this is actually 
quite remarkable, as both uiv± and U2V2 are separately comparable to on 
<9C 7 . Moreover, uiVi+U2V2\ x =o = 0(e~^ w / £ £~ 2 ). With these computations 
in mind, we observe that the general solution to (8.18) is given by 

U = C\U\ + C 2 U2 + U\Vl + u 2 v 2 - 

By our estimates on u±vi + ^i^lec-y an d the definitions of u\ and 112, we see 
that we may take c\ = 0(£) and C2 = 0{£ 2 ). Thus we compute that 

u(0) = c lUl (0) + c 2 n 2 (0) + m(0)t;i(0) + u 2 (0)v 2 (0) 
= c 2 n 2 (0) + u 2 {0)v 2 {0) since m(0) = 

= 0{£) + 0{e-^'H' z ) 
= 0(£). 

We conclude that -2(A - 2)~ 1 ^ = u < n(0) = 0{£) and so 



I -2(A - 2r 1 ^- = I ' uds< 0{£) f ds 

= 0{£ 2 ^ 



Combining the estimate for this term of the Hessian with the estimate for 
the other term of the Hessian discussed in Step 2 concludes the proof of the 
lemma. 

□ 



9. The Thurston metric and the Weil-Petersson metric 

From the formula (1.2), we can easily derive the result [Wol86b] (see also 
[McM08] and [Bon88]) that the Thurston metric is a multiple of the Weil- 
Petersson metric. 

9.1. The Thurston Metric. We begin by recalling the Thurston metric. 
To define this, imagine a sequence {7^} of closed curves which are becoming 
equidistributed in the sense that if B is a ball in the unit tangent bundle 
T 1 S', and we lift 7„ to its representative in T 1 ^, then 

£(-f n nB) Volume(S) 
™ £{ ln ) ~ Volume^ 1 ^)' 
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Thurston noted that since, for such a sequence of curves, -> as 

n — ► oo, then He ^^""- > would tend to a symmetric quadratic form on T(S); 
by the convexity of the length function, this tensor would be positive semi- 
definite, hence a (pseudo)-metric. Wolpert showed that 

Theorem 9.1. [Wol86b]. The Thurston metric is a multiple of the Weil- 
Petersson metric. 

The goal of the present section is to give a proof of this result that proceeds 
from evaluating formula (1.2) on a sequence {7 n } of curves that are becoming 
equidistributed in T l S. 



9.2. First Variation. We begin by first showing that dl(^ n ) / £{^ n ) — ► as 
n — > oo. Recall from (3.7) that 

d f Re$ 

-Mguin,t) = / ds. 



Ji VJH 1 10,1, J I 

at J 7n(0 ) 5o 

Now, the curves 7„ have unit tangent vectors ^7n(s) which equidistribute 
themselves in T X S, and so 

d L{g t ,-iri,t) = --J—z ( ^^-ds 



'in 90 

Re<3? ds 



ln 90 (-{in) 

f[[ Re<5>(p,9)dvol T i S (p,9) 
~" JJJ 90 vo\(T^S) 

where we interpret the meaning of the notation as follows. In the discussion 
so far, we have written Re $ to denote the value of the expression Re tp, when 
the quadratic differential^ = ipdz 2 was written in coordinates z = x + iy 
with being tangent to the geodesic. Now the vector field lifts to 
the canonical vector fields in T 1 S' tangent to the geodesic flow. We let 
the expression Re<£(p, 9) denote the value of Re<& on the surface in terms 
of a coordinate z = x + iy in which the geodesic direction described by 
(jp, 6) £ T l S is in the coordinate direction J^. Of course, if we change surface 

coordinates so that zq = e~ ie z, then for <E> = (pgdzg, we have ipg = e 2l0 (fo. 
Thus, when we integrate along the fiber of T 1 ^ — > S (with respect to 9 in 
the coordinates (z, 9) for T l S), we find / Re ^P^ d9 = 0. 
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9.3. The second variation. We recall the (second) formula (1.2) for the 
second variation of length: 

_i2|$ |2 



^/- <A - 2) " 



., ds 
% 

+ JT^ Im n(p) H > 2 / Im fl ( yq )ds(p)ds(q) 

ihn)JJy n xy n 2sinh(^) 



cosh(eZ(p, q) 
— 

|$|2 



-cis 



^/- 2(A - 2) ~^ 

1 /" Im$(p) f cosh{d{p,q) -%»)/2) Im$(g) 
t(ln)J, n 9o(p) J ln 2sinh(£( 7n )/2) 50 (g) 9P 

= A + / 2 - 

We examine the two integrals I\ and I2 separately. The first integral is 
immediate, as 

M 

by the equidistribution property, and 



M 



as the operator (A — 2) 1 is self adjoint. Then as — 2(A — 2) = 1, we 
find that 

, . ff |$| 2 dArea 

9.4. The Integral /2- We turn next to I2, where the computation is a bit 
more involved. Our basic plan mirrors our discussion of the first variation; 
we extend the terms in the integrand of I2 to all of T 1 S', and then integrate 
over the circular fiber to be left with an integral over the surface. 

We begin our more detailed discussion of the second (energy) integral by 
considering the version (4.33) of it in terms of a geometric kernel, i.e. 

1 f Im$(p) f cosh(d(p, q) — L/2) lm.$(q) 

v -dqdp 



Kin) J-yo 9o(p) J l0 2sinh(L/2) g (q) 
where L = £(70)- 

Now considering 70 as an embedded curve in the unit tangent bundle 
T 1 (»S, go), if we fix the point p = (p, v) G T 1 ^ as representing a point p G S 
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and a unit vector v 6 TpS, then a point g = (q,w) along 70 at distance t 
from p could be written 

q = exp ? to = G t p 

where Gt denotes the geodesic flow in T^S 1 for distance t. 
Thus we see that as L — > 00, the integral I2 converges to 

lim j _ f Im*(p) Z" 00 e-*^Im$(G t p) + Im$(G_ t p) A cftdp 



.Ms fob) 7o 2 9o(G t (p)) g (G- t (p)) vol^S)' 

Of course, as p varies in the fiber {(p, e t6 v)} over jj 6 S, we observe the 
points Gt(p,v) also arising as G-t(p, —v), and so we may rewrite the above 
limit integral as 

lim I — [ Im$ ^ [°° e -t l ™®( G tP) dtdp 
L ^oo 2 J Tls go (p) J go(G t (p)) vol^S)' 

To evaluate this last integral, imagine p = in the disk {\z\ < 1} and we 
represent the hyperbolic metric as 

5o=(*^=5o(0)(l-r 2 r 2 H 2 - 

An important matter here (as it was in the calculation of the first variation) 
is the question of how to interpret the meaning of Im$(g) in these coor- 
dinates: recall that we understood Im — to be the value of Im — when we 

go so 

defined the geodesic 70 as a vertical line in the coordinate system. In the 
fixed coordinate system of the disk {\z\ < 1}, write = 4>(z)dz 2 ; then since 
the hyperbolic geodesic through the origin and the point z = re tS is given 
by the ray te te , we see that Im$| re ie = Im e 2ie (f)(re td ). With this notation, 
our integral becomes 



lim I2 

L^oo 



r r2ir 

Js Jo 



Ime 2i V(0) e _t lme 2i9 <f>(r{t)e ie ) dtd6d Area 



sJo 9o(p) 50 (p)(l-r(t) 2 )- 2 27rArea( < S) 

f f°° e _ *(l -r(t) 2 ) 2 f 2n T /9?73l/ , NT , m ia „dtd8d Area, 

= / / 2, I / Im e 2, V Im e 2l ^ r t e' fl — — . 

JsJo 9o(P) Jo 2vrArea(S') 

Using the change-of-coordinates formula e~ % dt = 2(1 + r)~ 2 dr and the 
mean value theorem for harmonic functions (together with a half-angle for- 
mula to simplify the averaging), we find 



lim h= I — 1^ (T 2(1- rfdr 



1 0(0) 1 2 dArea 
2 Area(S) 
\$(p)\ 2 dAiea(p) 



Is 9o{p) 2 Area(S) ' 
Combining (9.1) and (9.2), we verify Theorem 9.1. In particular, 
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Hess£(7„) _ 4 r |<£>(p)| 2 dArea(p) 
n™ ro £( ln ) -3j s 50 (p)2 Area(S) 

4 2 
-3Area(5) llMllwp - 

Remark. The constant | found here agrees with that found by Wolpert 
[Wol86b] and McMullen [McM08]. See the comments ([McM08], p. 376) of 
McMullen on the consistency of conventions. 
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